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SECTION  1 
INTRODUCTION 

Although  the  literature  is  replete  with  analytical  studies  of  the  linear  dynamic  re¬ 
sponse  of  submerged  structures,  the  dynamic  instability  of  such  structures  has  received 
relatively  little  attention  [1],  In  1965,  Di  Maggio  [2]  studied  the  unstable  dynamic 
response  of  an  infinite  flat  plate  with  a  sinusoidal  imperfection  in  one  direction  sub¬ 
jected  to  an  in-plane  static  loading  in  that  same  direction;  the  plate  was  suddenly  re¬ 
leased  so  as  to  interact  with  an  acoustic  medium  on  one  side  of  the  plate.  He  found  that, 
in  the  vast  majority  of  cases,  the  acoustic  medium  may  be  treated  in  the  incompressible 
approximation.  In  1972,  Deng  and  Popelar  [3]  studied  the  parametric  instability  of  a  sub¬ 
merged  cylindrical  shell  initially  undergoing  sinusoidal  breathing  motions.  They  also 
found  that  the  acoustic  medium  could  be  accurately  treated  as  imcompressible .  References 
[4-6]  report  analyses  of  dynamically  excited,  submerged  shells  that  exhibit  instability 
characteristics.  In  all  of  these,  however,  approximate  treatments  of  the  fluid-structure 
interaction  are  used,  which  raises  questions  regarding  the  accuracy  of  the  results. 

This  report  presents  a  rigorous  treatment  of  the  dynamic  instability  of  an  infinite, 
elastic,  circular  cylindrical  shell  excited  by  a  transverse,  transient  acoustic  wave.  The 
field  quantities  appearing  in  appropriate  kinetic-energy,  potential-energy,  and  work-poten¬ 
tial  functionals  are  expanded  in  circumferential  Fourier  series, and  high-order  terms  are 
eliminated  in  a  consistent  manner.  The  residual  potential  formulation  [7-9],  which  consti¬ 
tutes  an  exact  formulation,  i s  used  to  treat  the  fluid-structure  interaction.  The  result¬ 
ing  modal  response  equations  provide  a  complete  and  rigorous  description  of  the  dynamic 
processes . 

The  modal  response  equations  are  integrated  numerically  in  time  for  excitation  by  plane 
waves  of  rectangular  and  exponential  pressure-profile.  Transient  response  histories  are 
provided  that  portray:  1)  modal  response  as  a  function  of  incident-wave  magnitude,  2)  the 
effects  of  ambient  hydrostatic  pressure  and  "live-load"  forcing  terms  (i.e.,  terms  that 
account  for  finite  translations  and  rotations  of  the  shell),  3)  shell  response  at  various 
locations,  4)  the  impact  of  flexural  stiffness  on  mode  participation,  and  5)  the  accuracy 
of  the  doubly  asymptotic  approximation  [8,10]  for  treatment  of  the  fluid-structure  interac¬ 
tion  in  dynamic  elastic  instability  analysis.  Conclusions  drawn  from  the  examination  of 
these  histories  are  listed  in  Section  TV. 
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The  present  report  may  be  regarded  as  contributing  to  the  second  stage  of  a  three- 
stage  study  of  the  transient  excitation  of  submerged,  infinite  cylindrical  shells.  The 
first  stage,  exemplified  by  References  [7]  and  [11],  deals  with  linear-elastic  shell  re¬ 
sponse.  The  second  stage,  exemplified  by  this  report,  deals  with  nonlinear-elastic  shell 
response.  The  last  stage,  exemplified  by  Reference  [12],  deals  with  nonlinear-inealstic 
shell  response.  A  corresponding  study  for  spherical  shells  would  be  most  useful,  the  first 
stage  of  which  is  exemplified  by  References  [13]  and  [10]. 


SECTION  2 

GOVERNING  EQUATIONS 


Consider  the  two-dimensional,  plane-strain  motions  of  the  submerged,  infinite,  cir¬ 
cular  cylindrical  shell  shown  in  Figure  1.  The  shell  is  excited  by  a  transient  acoustic 
wave  that  first  contacts  the  shell  at  8  =  it:  The  shell  is  thin  and  remains  elastic  at  all 
times;  geometric  nonlinearity  is  considered,  however,  which  introduces  the  possibility  of 
dynamic  instability. 

2.1  ENERGY  EXPRESSIONS 

Kinetic  and  strain  energy  expressions  for  the  cylindrical  shell  of  Figure  1  are: 

2tt 

T  =  pQh  J  (v2  +  w2)  a  dp 
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2tt  h/2 
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where  v  =  dv/dt,  etc.,  is  the  circumferential  stress,  e  is  the  circumferential 

0  0 

strain  and  z  is  the  thickness  coordinate.  Stress-strain  and  strain-displacement  rela¬ 
tions  for  the  shell  are  [14] 
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where  only  first-order  geometric  nonlinearities  are  considered.  The  introduction  of  (2) 
into  (1)  then  yields 
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An  expression  for  the  work  potential  appropriate  to  a  pressure  field  acting  on 
the  surface  of  a  smooth  shell  has  been  provided  by  Cohen  [15].  Fur  the  present  problem, 
that  expression  reduces  to 


2tt 

IT  =  /  [p(a  +  -^  +  w)  w  +  |-  (a  -  p)  w2  +  |-pv2  +  vw^]  do  (A) 

o  r  3 

where  p  =  p(r,0,t)  denotes  the  total  pressure  field.  Although  Cohen  lists  continuity 
of  the  pressure  field  as  a  requirement  for  the  existence  of  the  work  potential,  it  is 
easily  shown  that  a  more  lenient  requirement  is  satisfactory,  viz.,  that  the  pressure 
field  contain  a  finite  number  of  integrable  discontinuities. 

Now  the  displacement  and  pressure  fields  may  be  expanded  in  Fourier  series  as 
follows : 
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Also,  the  n  ^  0  Fourier  coefficients  for  v  and  w  may  be  transformed  into  extensional 
and  flexural  coefficients  as  follows  [9]: 
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The  incorporation  of  (5)  and  (6)  into  (1),  (3)  and  (4)  then  yields 
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where  p*  =  p  (a,t),  9pa/3r  =  [dp  (r,t)/3rj  ,  g  =  2ne  +  [(n2-l)/n]  f  ,  e  =  p  =  0 

n  n  n  n  r-a  n  n  n  m-n  m— n 

for  m  <  n,  and  the  6 ^  are  Kronecker  deltas.  Equations  7  constitute  the  basic  energy 
expressions  required  for  the  present  study. 

2.2  ELIMINATION  OF  HIGH-ORDER  TERMS 


It  is  now  appropriate  to- eliminate  from  (7)  those  terms 
necessary  for  a  consistent  formulation.  For  this  purpose,  the 
taken  to  be  of  the  order  of  the  critical  buckling  pressure  for 


of  order  higher  than  that 
n  =  0  pressure  harmonic  is 
the  shell,  which  is  [16] 
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Hence  pQ  ~  E(h/a)  ,  so  that  the  two  singly  underlined  terms  in  (7),  which  govern  linear, 

2 

static,  axisymmetric  response,  yield  wQ/a  ~  (h/a)  .  Next,  flexural  displacements  are  taken 

as  f  /a  ~  h/a,  so  that  the  doubly  underlined  terms,  which  govern  linear,  static,  flexural 
n  a  4 

response,  yield  p  ~  E(h/a)  .  Finally,  the  triply  underlined  terms,  which  govern  linear, 

n  3 

static,  nonaxisymmetric-extensional  response,  yield  en/a  ~  (h/a)  .  The  use  of  these  order- 

or-magnitude  relations  in  (7),  followed  by  the  elimination  of  terms  of  order  (h/a)^  and 

higher,  yields  the  simplified  energy  expressions 
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where  a  =  (n  -l)/n  and  A,  „  =  n,  a.a  a  [5  6  ,  .  +  6  „  ,  .  -  6  ,,  ,  .  .  + 

n  kimn  k  l  m  n  (k-S.)o  (m-n)o  (k-i.)  (m-n)  (k-Jt)(mfn) 

^(k+£) (nrt-n) ^ *  Note  that,  because  =  0,  the  lower  limits  of  six  of  the  summations  in  (9) 
may  be  changed  from  1  to  2. 

A  few  remarks  about  (9)  are  in  order.  First,  nonaxisymmetric-extensional  response 
fails  to  appear  as  a  significant  energy  contributor.  Second, in  the  reduction  of  (7)  to 
(9),  it  has  been  assumed  that  a|3pa/3r|  ~  |p^|  for  all  n.  Third,  the  last  of  (9)  contains 
no  nonlinear  terms  involving  pn;  this  implies  that  the  last  term  in  the  integrand  of  (4)  is 
unimportant.  Fourth,  a  nondimensional  parameter  y  has  been  introduced  into  the  flexural 
strain-energy  term  to  permit  consideration  of  a  sandwich  shell  consisting  of  two  concentric 
shells  of  thickness  h/2  separated  by  a  uniform  core  of  negligible  mass  and  in-plane  stiff¬ 
ness.  Such  a  shell  serves  as  a  convenient  plane-strain  model  for  a  stiffened  shell  [7,  17]. 

2  2  1 

With  D  as  the  sandwich  shell's  flexural  stiffness,  y  =  12(l-v  )D/Eh  ,  Finally,  the 
simplified  expression  for  the  shell's  strain  energy  is  positive-definite;  this  is  readily 
seen  by  observing  that,  with  w(e,t)  *  wQ(t)  +  f(8,t),  it  constitutes  a  Fourier-series  de¬ 
composition  of  the  expression  [cf.  (3)] 
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de 
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where  a  prime  denotes  a  6-derivative  and  the  asterisk  a  6-integral.  This  corresponds  to 
the  use  of  Rayleigh's  inextensibility  assumption,  i.e.,  w  »  -3v/36  for  nonaxisymmetric 
shell  response  [18]. 

It  is  interesting  to  examine  results  produced  by  (9)  in  certain  special  circum¬ 
stances.  First,  consider  the  nonaxisymmetric,  linear,  free  vibrations  of  a  hydrostatic¬ 
ally  pressurized  shell.  In  this  case,  p  (r,t)  ■  P„,  p  (r,t)  ■  0  for  n  t  0,  and  the 
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flexural  displacements  are  infinitesimal.  The  application  of  Lagrange's  equation  [19] 
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where  L  •  T  -  U  -  JI,  then  yields,  for  q  -  wq, 
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For  q  =  f  ,  the  application  of  Lagrange's  equation  yields 
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p 
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(13) 


The  Introduction  of  (12)  into  (13),  followed  by  the  assumption  of  sinusoidal  free  vibra¬ 
tion,  yields  the  modal  natural  frequency  equation 


<l-VPCn> 


(14) 


where  c  =  E/p  (1-v  )  is  the  plate  velocity  for  the  shell  material  and  =  p  c  (h/a)- 
2  o  2  °  2  o  o 

(n  -1)  (yh)  /12a  is  the  critical  pressure  for  the  nth  flexural  mode;  note  that  P^  = 

[cf.  (8)].  Equation  (14)  clearly  corresponds  to  the  flexural  frequency  equation  for  a 

pressurized  ring  [ 20] • 


Next,  consider  the  response  of  an  unpressurized  shell  to  nearly  uniform  radial 
impulse-excitation.  The  application  of  Lagrange's  equation,  (11),  to  the  simplified 
energy  expressions,  (9) ,  yields  for  axisymmetric  extensional  and  nonaxisymmetric  flex¬ 
ural  response 
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where  A,^  =  A,  +  A,  +  A,  +  A  and  the  discussion  following  (9)  has  been 

kirn  kiran  kinm  knim  nkira 

utilized. 

These  equations  are  similar  to  equations  presented  in  [21-23],  which  treat  this  partic¬ 
ular  problem  in  considerable  detail.  Equations  (15)  are  superior,  in  fact,  to  the  corres¬ 
ponding  equations  in  those  references  because  only  (15)  exhibit  all  of  the  following 
characteristics:  1)  n=l  rigid-body  motion  is  decoupled  from  n«0  breathing  motion. 
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2)  linear-vibration  frequencies  for  the  flexural  modes  are  given  by  the  first  of  (14)  with 

P  =0,  and  3)  the  associated  strain  energy  expression  is  positive-definite. 

H 

From  this  brief  examination  of  two  special  cases,  it  is  concluded  that  (9)  are 
suitable  for  the  present  analysis.  This  suitability  derives  from  the  consistent  elimina¬ 
tion  of  high-order  terms  from  appropriate  energy  expressions. 

2.3  FLU ID -STRUCTURE  INTERACTION 

A  rigorous  formulation  of  the  fluid-structure  interaction,  which  must  be  considered 
in  conjunction  with  (9),  may  be  constructed  as  follows  [7-9].  First,  fluid  pressure  and 
radial  fluid-particle  velocity  are  expressed  as  derivatives  of  a  fluid  velocity  potential 
as 


P  =  -  P  $ 


u  =  dq>/dr 


(16) 


Second,  the  total  acoustic  field  is  treated  as  the  superposition  of  an  acoustic  field  for 
the  (known)  incident  wave  and  an  acoustic  field  for  the  (unknown)  scattered  wave,  i.e., 

<p(r,G,t)  *  (r,  0  ,  t)  +  <ps  (r,  0,  t)  (17) 

Third,  compatibility  of  radial  fluid-particle  velocity  and  radial  shell  velocity^  is  en¬ 
forced  at  the  wet  surface  of  the  shell  as 


u  (a,  p  ,  t)  =  w  (0  ,  t)  (18) 

Finally,  the  wave  equation  and  radiation  condition  for  each  circumferential  harmonic  of 
the  scattered  wave  [see  (5)]  are  replaced  by  the  equivalent  residual-potential  relation 
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in  which  the  residual  potential  sp^  is  given  by  the  convolution  relation 

t 

(pj^  (r,  t)  =  -  f  rn  (r,  t')  <pgn  (r,  t-t')  dt  (20) 

o 

where  the  r^  are  characteristic  functions  that  resemble  step-exponential  functions  [7]. 

Equations  16-19  may  now  be  utilized  to  produce,  for  each  circumferential  harmonic, 
the  flu id -structure  interaction  relations 
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where  tp^  =  tp^fa,  fc)  is  obtained  from  (20).  These  equations  constitute  the  optimum 
form  of  the  information  required  for  a  rigorous  treatment  of  the  fluid-structure  inter¬ 
action. 


2. A  MODAL  RESPONSE  EQUATIONS 

Convenient  nondimens ional  equations  may  be  obtained  through  introduction  of  the 
following  convention: 


w  =  w/a  , 


t  =  ct/a 


Application  of  this  convention  to  (9)  yields 
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where  w  =  dw  /dt,  f  =  df  /dt,  p'  =  [3p  /3rl~  ,  ,  and  the  discussion  following  (9)  has  been 
o  o  ’  n  no  ro  Jr=l 

utilized.  In  like  fashion,  nondimensionalization  of  (21)  and  (20)  yields 
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where  the  asterisk  denotes  temporal  convolution.  Note  that  all  nondimensional  acoustic- 
field  quantities  in  (23)  and  (24)  pertain  to  the  wet  surface  of  the  shell. 

The  application  of  Lagrange's  equation,  (11),  to  (23),  followed  by  the  appropriate 
utilization  of  (24)  and  the  second  of  (6)  with  e^  neglected,  yields  the  nondimensional 
modal  response  equations  for  the  submerged  shell 
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where  all  circumflexes  have  been  dropped,  and  where  p  ■*  PQh/pa,  8  **  (CQ/C)  an^  ^  “  (Yh)  / 
12a2.  The  nondimensional  critical  buckling  pressure  for  the  shell  is,  from  (8),  Pc  * 

3  u  6 

Equations  (25)  constitute  the  modal  response  equations  needed  for  a  rigorous  analy¬ 
sis  of  the  nonlinear  dynamic  response  of  a  submerged,  infinite,  elastic,  circular,  cylin¬ 
drical  sandwich  shell  excited  by  a  transverse,  transient  acoustic  wave.  They  lend 
themselves  to  step-by-step  numerical  integration  in  time,  producing  modal  response  his¬ 
tories  that,  through  (6)  with  e^  neglected  and  through  (5),  yield  corresponding  shell  res¬ 
ponse  histories.  From  (2),  (5),  and  (6)  with  neglected,  extensional  and  flexural 
strain  response  histories  may  be  obtained  as 
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where  e  is  nondimensional  [see  (22)].  For  the  sandwich  shell  described  after  (9),  the  dis¬ 
tance  from  the  neutral  axis  to  the  outer  shell  fiber  is  given  by 
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SECTION  3 
NUMERICAL  RESULTS 

The  numerical  results  presented  in  this  section  have  been  generated  by  the  application 
of  fourth-order  Runge-Kutta  numerical  integration  to  the  first  four  of  (25)  and  the  use  of 
trapezoidal  integration  for  the  last  of  (25) .  A  variable  incrementing  procedure  with 
0. 002< At  <  0.2,  has  been  used.  Results  have  been  obtained  for  two  steel  shells  submerged  in 
sea  water,  all  characterized  by  PQ/p  =  7.72,  cq/c  =  3.53  and  a/h  =  100.  The  two  shells 
differ  in  terms  of  their  y-values,  which  are  5  and  10;  these  values  correspond  to  moderate 
and  heavy  stiffening  of  a  uniform  homogeneous  shell.  The  shells  are  excited  by  plane  acous¬ 
tic  waves  of  rectangular  and  exponential  pressure  profile  that  make  initial  contact  at  time 
t  =  0  along  the  line  0  =  n.  The  generalized  excitation  functions  for  these  waves  are  given 
by 

IT 

<PIn(t)  =  (_l)n+1  -2.  pj  y*  g(t  -  1  +  cos  C)  cos  nc  dc 
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(28) 
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uT  (t)  =  (-l)n+1  —  PT  f  g(t-l+cos?)  cose  cosne  d e 

in  tt  i  J 

o 

where  e  =1  for  n  =  0  and  e  =2  for  n  >  1,  and  PT  is  the  maximum  value  of  the  incident- 
n  n  ’I 

wave  pressure  profile.  For  the  rectangular  wave,  g(t)  =  H(t)  -  H(t-T),  where  H(t)  is  the 
Heaviside  step-function  and  T  is  the  pulse  duration;  for  the  exponential  wave,  g(t)  = 

H(t)  e-*^,  where  \  is  the  decay  constant.  Note  that  all  of  the  preceding  quantities  are 
nondimens ional,  having  been  normalized  in  accordance  with  (22). 

3.1  MODAL  RESPONSE 

Figure  2  shows  displacement  response  histories  for  the  (n«0)  breathing  mode  of  the 
y  =  10  sandwich  shell  when  excited  by  broad,  rectangular,  incident  waves  of  duration  T  - 
10.  The  pressure  magnitudes  of  these  waves  vary  from  1%  to  300%  of  the  critical  buckling 
pressure  for  the  shell  [see  (8),  or  the  discussion  following  (25)].  Note  that  all  re¬ 
sponses  are  normalized  to  the  magnitude  of  the  incident  wave,  so  that  coalescence  of  re¬ 
sponse  histories  implies  linearity  of  response.  In  this  connection,  it  is  seen  that  the 
response  for  P^  -  P^,  and  =  2P^  are  virtually  coincident  with  the  (linear)  response  for 
P  »  0.01  P  •  nonlinear  effects  are  discernible  for  PT  *  3P„,  but  they  are  relatively  un- 
important. 
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Velocity  response  histories  for  the  (n=l)  translational  mode  are  shown  in  Figure  3. 
Nonlinear  effects  are  barely  discernible.  This  implies  that  the  last  term  in  the  last  of 
(23)  is  of  negligible  importance  because,  without  it,  the  n=l  mode  is  totally  uncoupled 
from  the  other  modes  and  is  governed  by  purely  linear  equations  [see  (25)  and  the  second 
of  (24)]. 


The  n=2  lobar  mode  is  the  one  that  exhibits  significant  nonlinear  behavior,  as  shown 
by  the  displacement  response  histories  of  Figure  4.  The  response  may  be  conveniently  de¬ 
scribed  as  occurring  in  four  phases:  an  "envelopment  phase",  which  extends  from  t=0  to 
t=2,  a  "pressurization  phase",  which  extends  from  t=2  to  t=T,  a  "development  phase", 
which  extends  from  t=T  to  t=T+2,  and  a  "f ree-vibration  phase",  which  extends  onward  from 
t=T+2.  During  the  envelopment  phase,  in  which  the  incident  wave  front  is  passing  over  the 
shell,  the  response  is  essentially  linear.  The  pressurization  phase,  in  which  the  shell  is 
essentially  hydrostatically  pressurized  by  the  incident  wave,  is  characterized  by  either 
oscillatory  or  exponential  response,  depending  upon  the  magnitude  of  Pj,  During  the  devel¬ 
opment  phase,  in  which  the  back  of  the  rectangular  incident  wave  is  passing  over  the  shell, 
the  response  exhibits  sudden,  but  modest  reversal.  Finally,  the  free-vibration  phase  is 
characterized  by  low-frequency  sinusoidal  motion.  Clearly,  the  appearance  of  response 
overshoot  during  the  free  vibration  phase  depends  upon  exponential  growth  experienced  dur¬ 
ing  the  pressurization  phase,  which,  in  turn,  depends  upon  the  magnitude  of  P^. 

Displacement  response  histories  for  the  n=3  lobar  mode  are  shown  in  Figure  5.  The 
preceding  description  of  n=2  response  is  applicable  here  also.  Especially  visible  in  the 
n=3  response  histories  is  the  virtually  undamped  nature  of  the  sinusoidal  motions  during 
the  free-vibration  phase.  This  indicates  that  the  surrounding  fluid  provides  very  small 
acoustic  radiation  damping  for  these  motions,  which  is  to  be  expected  when  the  characteris¬ 
tic  structural  wavelength  2ira/n  is  much  smaller  than  the  acoustic  wavelength  c/f,  where  f 
is  the  frequency  of  oscillation  [1],  The  n=4  and  n=5  lobar  modes  have  been  included  in 
these  y  =  10  shell  calculations,  but  exhibit  peak  displacements  substantially  smaller  than 
that  of  the  n=3  mode.  Hence  response  histories  are  not  shown  for  these  modes. 


It  is  instructive  to  examine  dynamic  instability  of  the  flexural  modes  during  the 
pressurization  phase.  During  this  phase,  n=0  displacement  is  approximately  equal  to  the 
hydrostatic  value 
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Also,  flexural  response  is  relatively  slow,  so  that  <Pgn  «  <pgn  and,  from  the  last  of  (25) 
and  [7], 
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Hence  the  fourth  of  (25)  yields,  with  u^n  =  0  for  n  >  2,  2  <  t  <  T  [9], 


j  n 


l 

n 


f 

n 


(30) 


(31) 


Equations  (29)  and  (31),  along  with  the  approximations  £  =»  0,  wq  «  0  and  neglect 

of  the  triple  summation  in  the  second  of  (25),  yield  for  flexural  response  during  the 
pressurization  phase 


<u  S-P  +  ;>"„  +  ("2  -  1)  <PC„  -  Pj)  f„  ■  0  (32) 

n 

2 

where  PCn  =  (n  -  l)pg£.  This  equation  exhibits  the  effect  of  added  fluid  mass  associated 
with  low-frequency  shell  response;  it  also  indicates  that  shell  response  is  oscillatory  for 
Pj  <  PCn  and  exponential  for  P^.  >  P  ,  as  seen  in  Figures  4  and  5.  It  is  worth  noting  that 
with  the  removal  of  P  ,  (32)  also  governs  flexural  mode  response  during  the  free-vibration 
phase. 


Modal  response  histories  for  excitation  by  incident  step-exponential  waves  with  decay 
constant  X  =  1  are  shown  in  Figures  6-9.  The  normalized  responses  are  substantially  smal¬ 
ler  than  their  rectangular-wave-excitation  counterparts  in  Figures  2-5  because  of  the  re¬ 
duction  in  pulse  width.  Otherwise,  the  behavior  is  similar,  with  significant  nonlinear 
effects  appearing  only  in  the  n  5  2  harmonics. 

3.2  AMBIENT  HYDROSTATIC  PRESSURE  AND  LIVE  LOAD 


It  is  interesting  to  examine  the  effects  of  ambient  hydrostatic  pressure  on  the  shock 
response  of  the  y  =  10  shell.  The  modifications  in  (25)  required  for  such  an  examina- 

o  o 

tion  merely  involve  the  replacement  of  wq  by  wq  -  P^/pg  and  of  by  <p^  -  PH,  where  PH  is 

the  magnitude  of  the  hydrostatic  pressure.  With  these  replacements,  the  first  of  (25)  re¬ 
mains  the  same,  except  that  it  now  pertains  to  breathing  motions  about  a  static  equilibrium 

2  2 

radius  of  1  -  PH/pg.  The  second  of  (25)  changes  only  to  the  extent  that  the  term  £(n  -1) 

is  now  multiplied  by  (1  -  Pu/P_  ),  where  P_  is  given  after  (32). 

n  on  on 

For  static  stability,  P„  must  be  less  than  the  smallest  of  the  P  ,  which  is  P„„  »  P„ 

H  on  o i  o 

*=  3pB£.  Hence  the  reduction  in  static  equilibrium  radius  cannot  exceed  3£,  which  is  very 
small.  Thus  the  principal  effect  of  ambient  hydrostatic  pressure  is  the  reduction  of  flex¬ 
ural  stiffness,  which  is  a  destabilizing  influence. 

It  is  also  interesting  to  examine  inaccuracies  introduced  into  transient  response  com¬ 
putations  by  the  neglect  of  the  "live- load"  terms  in  (4).  These  terms  account  for  the 
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effects  of  finite  translations  and  rotations  of  the  shell  on  the  work  done  by  the  normal 
pressure  loading  as  referenced  to  the  undeformed  shell  surface.  All  the  terms  in  (4)  ex¬ 
cept  p*a*w  are  live- load  terms;  when  processed  through  the  Fourier-decomposition  and  term- 
elimination  operations  of  Subsections  2.1  and  2.2,  they  appear  in  nondimensional  form  as 
the  last  two  summation  terms  in  the  last  of  (23).  Following  the  application  of  Lagrange's 

equation  (11)  and  the  introduction  of  the  fluid-structure  interaction  equations  (24),  live- 
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load  effects  manifest  themselves  as  the  terms  [  (n  -l)/n  ]  (<p ^  4-  <p  )  and  -wq  in  the  second 
of  (25). 

Response  computations  designed  to  demonstrate  the  effects  of  ambient  hydrostatic  pres¬ 
sure  and  live  load  have  been  performed  for  the  -y  =  10  sandwich  shell  excited  by  both  T  =  10, 
Pj  =  2?^  rectangular  and  X  =  1,  =  5P^  exponential  incident  waves.  As  would  be  expected, 

n=0  and  n=l  response  is  unaffected  by  the  introduction  of  ambient  hydrostatic  pressure  or 
the  omission  of  live-load  terms  in  (25).  In  contrast,  the  response  behavior  of  the  n=2 
flexural  mode  is  substantially  affected,  as  indicated  in  Figures  10  and  11.  It  is  clear 
from  these  figures  that  ambient  hydrostatic  pressure  and  live  load  are  both  significant  de¬ 
stabilizing  influences  for  this  mode*.  The  higher  modes  are  much  less  influenced,  however, 
as  indicated  by  the  n=3  response  histories  of  Figures  12  and  13. 

3.3  SHELL  RESPONSE 

Modal  response  histories  for  0  <  n  i  5  have  been  superposed  in  accordance  with  (5)  and 
(26)  to  construct  the  shell  response  histories  of  Figures  14-19.  Figures  14,  16  and  18 
pertain  to  excitation  by  a  T  =  10,  P^  =  2PC  rectangular  incident  wave,  while  Figures  15,  17 
and  19  pertain  to  excitation  by  a  X  =  1,  P^  =  5P^,  exponential  incident  wave. 

Shown  in  Figures  14  and  15  are  deformational  displacement  histories,  which  constitute 
displacement  histories  with  rigid-body  motion  removed,  i.e.,  w^(0,t)  =  w(0,t)  -  f^(t)  cos  6. 
It  is  seen  that  the  effects  of  ambient  hydrostatic  pressure  and  live  load  are  significant 
for  excitation  by  the  rather  broad  rectangular  pulse,  but  are  relatively  minor  for  excita¬ 
tion  by  the  rather  narrow  exponential  pulse. 

Figures  16  and  17  show  radial  velocity  histories  at  two  points  on  the  shell.  The  ef¬ 
fects  of  the  ambient  hydrostatic  pressure  and  live  load  are  clearly  negligible.  Strain  re¬ 
sponse  histories  are  shown  in  Figures  18  and  19  at  locations  selected  to  emphasize  flexural 
and  nonlinear-extensional  contributions  to  total  strain.  The  coincidence  of  the  w  -  and 
tg-histories  in  the  figures  demonstrates  that  the  nonlinear-extensional  term  in  the  first 


*  A  development  like  that  which  produced  (32),  but  which  includes  ambient  hydrostatic  pres¬ 
sure  and  excludes  live-load  terms,  leads  to  a  prediction  of  the  coalescence,  for  0<t^l0, 
of  the  response  histories  of  Figure  10  with  the  circle  and  triangle  designators. 
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of  (26)  is  minuscule.  This,  along  with  the  smallness  of  the  flexural  strain  history,  demon¬ 
strates  that  the  strain  response  of  the  shell  is  dominated  by  the  n=0  breathing  mode.  Hence, 
peak  strain  in  the  shell  is  accurately  estimated  as 

| e I  «  (P'1  I w  |  )  •  NP„  (33) 

'  1  max  I  1  o 1  '  C 

max 

where  N  =  PT/P  .  For  the  shell  of  Figure  2,  for  example,  P-^lw  |  =  1.125  and  P_  =  3v6C 

i  c  i  o  max  v> 

=  2.405  x  10“3;  hence  |s|max  ~  0.27N%,  as  suggested  by  Figure  18. 

3.4  FLEXURAL  STIFFNESS  EFFECTS 

It  is  informative  to  compare  the  response  behavior  of  a  moderately  stiffened  (y  =  5) 

shell  with  that  of  the  heavily  stiffened  (y  =  10)  shell  of  Figures  2-19.  As  the  flexural 

2 

stiffness  of  a  sandwich  shell  is  proportional  to  y  [see  the  discussion  preceding  (10)], 
the  P^-value  for  the  moderately  stiffened  shell  is  only  one-quarter  of  that  for  the  heav¬ 
ily  stiffened  shell.  Hence,  in  order  to  maintain  proper  calibration  between  the  excita¬ 
tion  levels  for  the  two  shells,  P^-values  of  0.04  Pc,  4  ?c,  8  Pc,  and  12  ?c  are  used  for 
the  moderately  stiffened  shell. 

Figures  20-24  show  modal  response  histories  for  the  y  =  5  shell.  A  comparison  of  the 
first  of  these  with  Figure  2  indicates  that  nonlinear  effects  in  n=0  response  are  more  pro¬ 
nounced  for  the  y  =  5  shell  than  they  are  for  the  y  =  10  shell;  even  so,  they  remain  rela¬ 
tively  unimportant.  Figure  21,  which  pertains  to  n=l  response,  is  virtually  a  duplicate  of 
Figure  3,  in  which  nonlinear  effects  are  barely  discernible. 

Figure  22,  which  pertains  to  n=2  response,  displays  significant  nonlinear  behavior,  as 
does  Figure  4  for  the  •,  =  10  shell.  For  the  same  value  of  P^.  >  Pc,  an  n=2  response  peak 
during  the  f ree-vibration  phase  for  the  y  =  5  shell  considerably  exceeds  its  counterpart 
for  the  y  =  10  shell.  This  is  suggested  by  (32),  which  predicts  exponential  response 
growth  for  P  >  P„  =  P  during  the  pressurization  phase.  For  example,  for  the  y  =  5 

^  2 

shell  with  =  12  P^,  the  stiffness  coefficient  multiplying  f ^  is  -3*11-PC  =  -99*yB*25h  / 

12a^,  which  exceeds  in  magnitude  its  counterpart  for  the  y  =  10  shell  with  P  =  3  Pr,  which 

2  2  -  it 

is  -3‘2‘P^  =  -18'u6"100h  /12a  .  In  contrast  to  the  P^  >  P  results  of  Figure  22,  the  peak 
f ree-vibration-phase  response  of  the  y  =  5  shell  for  P^.  =  0.04  Pc  is  only  one-half  that  of 
the  y  =  10  shell  for  P  =  0.01  Pc.  This  may  be  explained  as  follows.  An  approximate  ex¬ 
pression  for  the  linear ,  post-envelopment,  n=2  response  of  a  step-wave-excited  shell  is 
given  in  [24],  This  expression,  in  turn,  leads  to  the  following  expression  for  n*'2,  linear, 
f ree-vibration-phase  response  appropriate  to  excitation  by  a  rectangular  incident  wave 
whose  duration  T  is  much  smaller  than  the  n=2  mode's  natural  period  of  oscillation  during 
the  f ree-vibration  phase: 
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(34) 


4/5  1 

f2^  2\i  +  4/~5  PI  Q2  slnQ2  ^  ~ 

2 

where  =  36p8£/ (5p+2) .  Hence,  inasmuch  as  £  for  the  y  =  10  shell  is  four  times  that  for 
the  y  -  5  shell,  peak  free-vibration  response  of  the  y  =  10  shell  for  <<  is  double 
that  of  the  y  -  5  shell,  all  other  parameters  being  equal. 

It  is  clear  from  the  preceding  paragraph  that  the  nonlinear  response  of  a  flexural 
mode  to  rectangular-wave  excitation  depends  critically  upon  the  "initial  conditions"  for 
the  pressurization  and  free-vibration  phases  and  the  value  of  the  pressurization-phase 

"stiffness"  parameter  in  (32).  With  PT  =  NP„,  (32)  predicts  exponential  growth  during  the 

1  ^  A 

pressurization  phase  for  those  modes  whose  modal  index  satisfies  the  inequality  n  <  (3N+1)“. 

Hence  a  large  value  of  N  implies  that  many  flexural  modes  may  contribute  signif icantly  to 
shell  response,  while  a  small  value  of  N  implies  that  only  the  lowest  flexural  modes  need 

be  considered.  This  is  illustrated  in  Figures  23  and  24,  which  show  modal  response  his¬ 

tories  for  n=3  and  4.  Although  the  n=3  response  peaks  for  the  y  =  5  shell  are  no  larger 
than  their  counterparts  for  the  y  =  10  shell  shown  in  Figure  5,  the  n=4  mode  of  the  y  =  5 

shell  responds  strongly,  especially  for  P^  =  12  Pc.  This  is  in  contrast  to  n=4  response 

for  the  y  =  10  shell,  which  is  so  small  that  it  is  not  even  included  in  the  discussion  of 
Subsection  3.1. 

The  increased  participation  of  the  higher  flexural  modes  in  the  response  of  the  y  =  5 
shell  suggests  that  flexural  strain  might  now  play  a  much  more  important  role  than  that 
portrayed  in  Figure  18  for  the  y  =  10  shell.  This  is  not  the  case,  however,  as  shown  in 
Figure  25,  because  the  increase  in  flexural  mode  response  is  essentially  negated  by  the 
decrease  in  the  distance  between  the  outer  fibers  of  the  shell.  There  is  a  discernible 
nonlinear-extensional  contribution  to  total  strain,  but  it  is  hardly  significant. 

The  preceding  comparison  of  y  =  5  and  y  =  10  shell  response  illustrates  the  dissimilar¬ 
ities  between  the  response  behavior  of  moderately  stiffened  and  heavily  stiffened  shells  to 
a  given  incident  wave.  The  contrast  between  an  unstiffened  shell  and  a  moderately  stiffened 
shell  is  even  greater.  For  example,  Pj  =  8  P^  for  y  =  5  corresponds  to  P^  =  200  Pc  for 
y  =  1.  Thus,  while  (32)  predicts  exponential  growth  during  the  pressurization  phase  for 
the  n  <  5  flexural  modes  of  the  y  =  5  shell,  it  predicts  such  growth  for  the  n  <  25  flex¬ 
ural  modes  of  the  y  =  1  shell.  Experimental  observations  of  such  short-structural-wave¬ 
length  instability  are  reported  in  [25,  26J. 

3.5  DOUBLY  ASYMPTOTIC  APPROXIMATION 

As  mentioned  in  Section  1,  an  aspect  of  particular  interest  in  the  present  study  has 
been  the  accuracy  of  the  doubly  asymptotic  approximation  (DAA)  [8,  10]  for  treatment  of  the 
fluid-structure  interaction.  In  the  present  context,  such  treatment  involves  replacement 
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of  the  coefficient  £  in  the  third  and  fourth  of  (25)  by  n,  omission  of  ^  in  those  equa¬ 
tions,  and  omission  of  the  last  of  (25)  entirely.  With  these  changes,  modal  response  his¬ 
tories  have  been  calculated  for  direct  comparison  with  the  exact  histories  of  Figures  2-9 
produced  by  the  residual  potential  formulation  (RPF) . 

Modal  response  histories  pertaining  to  rectangular -wave  excitation  of  the  y  =  10  shell 
are  shown  in  Figures  26-29.  Two  of  the  four  histories  in  each  figure  pertain  to  the  larg¬ 
est  incident-wave  magnitude  considered  in  Figures  2-5,  while  the  other  two  pertain  to  an 
incident-wave  magnitude  so  small  that  linear  shell  response  is  assured.  An  examination  of 
Figures  26-29  reveals  that  discrepancies  between  the  DAA  and  RPF  histories  for  P^.  =  3  Pc 
are  similar  in  nature  and  size  to  those  between  the  DAA  and  RPF  histories  for  P^.  =  0.01  Pc. 
This  indicates,  of  course,  that  DAA  accuracy  for  elastic  dynamic  instability  calculations 
is  fully  comparable  to  that  for  linear-elastic  response  calculations.  In  particular,  the 
tendency  of  the  DAA  to  overestimate  acoustic  radiation  damping  for  low-frequency  free 
vibration  does  not  artificially  suppress  unstable  shell  response  during  the  pressurization 
phase. 

A  comparison  of  DAA  and  RPF  response  histories  for  exponential-wave  excitation  of  the 
Y  =  10  shell  leads  to  the  conclusions  expressed  in  the  previous  paragraph.  Hence  exponen¬ 
tial-wave  counterparts  to  Figures  26-29  need  not  be  shown.  A  thorough  examination  of  the 
accuracy  of  the  DAA  for  linear-elastic  shell  response  is  provided  in  127]. 


SECTION  4 


CONCLUSION 


This  study  has  dealt  with  the  dynamic  instability  of  an  infinite,  elastic,  submerged, 
circular  cylindrical  shell  excited  by  a  transverse,  transient  acoustic  wave.  Circumferen¬ 
tial  Fourier  decomposition  of  the  field  quantities  appearing  in  appropriate  energy  function¬ 
als,  followed  by  consistent  elimination  of  high-order  terms,  has  led  to  rather  simple  shell 
response  equations  with  satisfactory  attributes.  The  fluid-structure  interaction  has  been 
treated  rigorously  in  accordance  with  the  residual  potential  formulation,  which  has  been 
used  successfully  in  a  number  of  previous  studies. 

Numerical  results  have  been  presented  in  the  form  of  117  transient  response  histories 
pertaining  to  excitation  by  plane  waves  of  rectangular  and  exponential  pressure-profile. 
Examination  of  these  results  has  led  to  the  following  conclusions: 

1.  Dynamic  instability  effects  are  significant  only  with  respect  to  flexural  shell  re¬ 
sponse. 

2.  The  dynamically  unstable  flexural  response  that  occurs  while  the  incident-wave  profile 
passes  over  the  shell  profoundly  affects  subsequent  free-vibration  flexural  response. 

3.  For  peak  total  strains  smaller  than  1%,  only  the  lowest  flexural  modes  of  a  heavily 
stiffened  shell  experience  dynamic  instability;  if  the  degree  of  stiffening  is  reduced, 
however,  higher  flexural  modes  exhibit  such  instability  until,  in  the  limit  of  an  un¬ 
stiffened  shell,  numerous  flexural  modes  are  involved.  Hence,  the  response 

of  an  unstiffened  shell  to  a  given  incident  wave  is  significantly  different  than  that 
of  an  appreciably  stiffened  shell  to  the  same  wave. 

4.  n=0  breathing  motion  dominates  strain  response,  and  quasi-linear  (especially  n=l  trans¬ 
lational)  motions  dominate  velocity  response;  nonlinear  flexural  motions  are  important 
only  with  regard  to  deformational  displacement  response,  i.e. ,  displacement  response 
exclusive  of  rigid-body  translation. 

5.  As  a  consequence  of  the  preceding  conclusions,  live-load  forcing  terms  and  ambient 
hydrostatic  pressure  substantially  affect  deformational  displacement  response,  but 
have  minor  impact  on  velocity  and  strain  response. 

6.  The  live-load  terms  in  the  surface-pressure  work  potential  that  involve  pressure  grad¬ 
ients  may  be  neglected  [see  (4)]. 

7.  The  accuracy  of  the  doubly  asymptotic  approximation  for  treatment  of  the  fluid-struc¬ 
ture  interaction  in  dynamic  elastic-instability  calculations  is  fully  comparable  to 
that  observed  in  linear-elastic  response  calculations. 
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FIGURE  2.  n=0  Displacement  Response  of  the  y  =  10  Shell  to 

Rectangular  Incident  Waves  with  T  =  10. 
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FIGURE  8.  n=2  Displacement  Response  of  the  y  =  10  Shell  to 

Exponential  Incident  Waves  with  A  =  1. 
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FIGURE  11.  n=2  Displacement  Response  of  the  y  =  10  Shell  to  an 

Exponential  Incident  Wave  with  A  =  1,  P  =  5P~. 
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FIGURE  15.  Deformational  Displacement  Response  of  the  y  =  10  Shell  to 
an  Exponential  Incident  Wave  with  A  =  1,  PT  =  5P_. 
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FIGURE  16.  Radial  Velocity  Response  of  the  y  =  10  Shell  to  a 
Rectangular  Incident  Wave  with  T  =  10,  PT  =  2P„. 
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FIGURE  27.  Approximate  and  Exact  n“l  Velocity  Response  of  the  Y 
Shell  to  Rectangular  Incident  Waves  with  T  *  10. 


FIGURE  28.  Approximate  and  Exact  n*2  Displacement  Response  of 
the  y  -  10  Shell  to  Rectangular  Incident  Waves  with 
T  -  10. 


FIGURE  29.  Approximate  and  Exact  n*3  Displacement  Response  of 
the  y  ■  10  Shell  to  Rectangular  Incident  Waves  with 
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Naval 

ATTN: 

dama-css-n 

2  cy 

Harry  Diamond  Laboratories 

2  cy 

Department  of  the  Army 

2  cy 

ATTN: 

DELHD-l-TL 

ATTN: 

DELHD-N-P 

Naval 

U.S.  Army  Ballistic  Research  Labs 

ATTN: 

DRDAR-TSB-S 

U.S.  Army  Engr.  Waterways  Exper,  Station 

2  cy 

ATTN: 

Library 

Naval 

U.S.  Army  Material  &  Mechanics  Rsch.  Ctr, 

ATTN: 

DRXMR-TE,  R.  Shea 

Naval 

U.S.  Army  Nuclear  &  Chemical  Agency 

ATTN: 

Library 

Naval 

DEPARTMENT  OF  THE  NAVY 

ATTN 

ATTN 

ATTN 

ATTN 

ATTN 

ATTN 

ATTN 

ATTN 


Code  1740.1 
Code  173 
Code  1844 
Code  1740.6 
Code  L42-3 
Code  11 
Code  172 
Code  2740 
Code  1770.1 
Code  1740.5 


ATTN:  Code  09M22C 


ATTN:  MAT  08T-22 


ATTN:  Code  4471 


ATTN: 

ATTN: 


Code  0142 
Code  69NE 


Naval  Research  Laboratory 


ATTN 

ATTN 

ATTN 

ATTN 

ATTN 

ATTN 

ATTN 

ATTN 


Code  8003 
Code  8301 
Code  6380 
Code  8100 
Code  8406 
Code  8445 
Code  2627 
Code  8440 


Naval  Sea  Systems  Comnand 


ATTN 


SEA-09G53 

SEA-3221 

SEA-08 

SEA-323 


Naval  Surface  Weapons  Center 


ATTN 

ATTN 

ATTN 

ATTN 


Code  R13 
Code  F31 
Code  RIO 
Code  R15 
Code  R14 


ATTN:  Tech.  Library  &  Info.  Services  Branch 


ATTN:  Code  233 


ATTN:  Code  10 


Naval  Construction  Battalion  Center 
Civil  Engineering  Laboratory 
ATTN:  Code  L08A 

Naval  Electronic  Systems  Command 
ATTN:  PME  117-21 


New  London  Laboratory 

Naval  Underwater  Systems  Center 

ATTN:  Code  401,  J.  Kalinowski 
ATTN:  Code  401,  J.  Patel 
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DEPARTMENT  OF  THE  NAVY  (Continued) 

DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 

Newport  Laboratory 

General  Electric  Company-TEMPO 

Naval  Underwater  Systems  Center 

ATTN:  DAS I AC 

ATTN:  Code  EM 

ATTN:  Code  363,  P.  Paranz ino 

University  of  Illinois 

Office  of  Naval  Research 

Consulting  Services 

ATTN:  N.  Newmark 

ATTN:  Code  715 

2  cy  ATTN:  Code  474,  N.  Perrone 

Institute  for  Defense  Analyses 

ATTN:  Classified  Library 

Office  of  the  Chief  of  Naval  Operations 

ATTN:  OP  982 

Kaman  AviDyne 

ATTN:  OP  953 

ATTN:  Library 

ATTN:  OP  981N1 

ATTN:  OP  957E 

Kaman  Sciences  Corp. 

ATTN:  OP  951 

ATTN:  Library 

ATTN:  OP  37 

ATTN:  OP  604C 

Lockheed  Missiles  and  Space  Co.,  Inc. 

ATTN:  OP  03EG 

ATTN:  Technical  Information  Center 

ATTN:  OP  987 

15  cy  ATTN:  T.  Geers 

ATTN:  OP  21 

ATTN:  C.  Yen 

Strategic  Systems  Project  Office 

Merritt  CASES,  Inc, 

Department  of  the  Navy 

ATTN:  Library 

ATTN:  NSP-272 

ATTN:  NSP-43 

Pacifica  Technology 

2  cy  ATTN:  J.  Kent 

DEPARTMENT  OF  THE  AIR  FORCE 

Air  Force  Institute  of  Technology 

Patel  Inteprises,  Inc. 

ATTN:  M.  Patel 

Air  University 

ATTN:  Library 

Physics  Applications,  Inc. 

ATTN:  C,  Vincent 

Air  Force  Weapons  Laboratory,  AFSC 

ATTN:  SUL 

R  &  D  Associates 

DEPARTMENT  OF  ENERGY  CONTRACTORS 

ATTN:  Technical  Information  Center 
ATTN:  C.  MacDonald 

Lawrence  Livermore  Laboratory 

SRI  International 

ATTN:  Document  Control  for  Technical 

ATTN:  G.  Abrahamson 

Information  Dept.  Library 

ATTN:  A.  Florence 

Los  Alamos  Scientific  Laboratory 

Tetra  Tech,  Inc. 

ATTN:  Document  Control  for  D.  Nowlin 

ATTN:  Library 

ATTN:  Document  Control  for  MS  364 

ATTN:  L.  Hwang 

(Class.  Reports  Lib.) 

Sandia  Laboratories 

Weidlinger  Assoc. ,  Consulting  Engineers 

3  cy  ATTN:  M.  Baron 

ATTN:  Document  Control  for  3141 

Sandia  Laboratories 

Weidlinger  Assoc.,  Consulting  Engineers 

ATTN:  J,  Isenberg 

Livermore  Laboratory 

ATTN:  Document  Control  for  Library  & 

Security  Classification  Division 

DEPARTMENT  OF  DEFENSE  CONTRACTORS 


BDM  Corp. 

ATTN:  T.  Neighbors 

Bolt  Beranek  &  Newman,  Inc. 
ATTN:  R.  Haberman 


Cambridge  Acoustical  Assoc.,  Inc. 
ATTN:  M.  Junger 


General  Dynamics  Corp. 
Electric  Boat  Division 
3  cy  ATTN:  M.  Pakstys 
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